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Abstract 

o ■ 

s ! ' This article provides a concise review of the problem of neutrino radiation from dense matter. 

The subjects addressed include quantum kinetic equations for neutrino transport, collision integrals 
describing neutrino radiation through charged and neutral current interactions, radiation rates from 
. pair-correlated baryonic and color superconducting quark matter. 

a\ ; 

^ 1 Introduction 

■ After an initial phase of rapid (of the order of weeks to years) cooling from temperatures T ~ 50 MeV 
down to 0.1 MeV, a neutron star's core settles in a thermal quasi-equilibrium state which evolves slowly 
over the time scales 10 3 — 10 4 yr down to temperatures T ~ 0.01 MeV. The cooling rate of the star 
t-h ■ during this period is determined by the processes of neutrino emission from dense matter, whereby the 
= neutrinos, once produced, leave the star without further interactions. Understanding the cooling processes 
[ that take place during this neutrino radiation era is crucial for the interpretation of the data on surface 
"^h 1 temperatures of neutron stars. While the long term features of the thermal evolution of neutron stars 
O !~ are insensitive to the initial rapid cooling stage, the subsequent route in the temperature versus time 
q diagram, which includes the late time (t > 10 5 yr) photoemission era, strongly depends on the emissivity 
?-) | of matter during the neutrino cooling era. 

This lecture is a concise introduction to the physics of neutrino radiation from dense nucleonic and 
, quark matter in compact stars. It starts with a classification of the reactions in Sec. 11.11 which is followed 
. ^ | by a discussion of quantum kinetic equations for neutrinos and neutrino emissivities in Sec. El In Sec. [3] 
^ 1 examples are given of polarization tensors of superfluid nucleonic matter and color superconducting quark 
^ . matter. We close by suggesting two exercises for students. 

_ _ i 

1.1 Classification of the reactions 

Historically, the weak reactions in neutron stars were classified within the quasiparticle description for 
fermions in matter: each reaction is distinguished by the number of the participating quasiparticles and 
the weak-interaction current. The simplest neutrino emission processes that involve single fermionic 
quasiparticle in the initial (final) state can be written as 

fi - h + e + v, h + e - h + v, (1) 
/ -> / + v + v, (forbidden) (2) 



* Based in part on the lectures delivered at the Summer school "Dense Matter In Heavy Ion Collisions and Astro- 
physics", BLTP, Joint Institute for Nuclear Research, Dubna, Russia. The complete lectures are available online at 
http/ / theor . j inr . ru/ ~dm2006/talks . html 
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where the first line is the charged current /3-decay and its inverse, with fi and f2 being neutron and 
proton quasiparticles in nucleonic matter or d and u flavor quarks in deconfined quark matter; / refers 
to a fermion. This process is known in astrophysics as the Urea processes pp. The Urea reaction is 
kinematically allowed in nucleonic matter under /3-equilibrium if the proton fraction is sufficiently large, 
Y p > 11— 14% [2j. In deconfined, chirally symmetric, and interacting quark matter at moderate densities 
the Urea processes is kinematically allowed for any asymmetry between u and d quarks [3]. The second 
process - the neutral current neutrino pair bremsstrahlung, Eq. ([2]), is forbidden by the energy and 
momentum conservation, if one adopts the quasiparticle picture. If, however, we choose to work with 
excitations that are characterized by finite widths, the reaction ([2]) is allowed [I]. The processes with two 
fermions in the initial (and final) states are the modified Urea and its inverse [5] 

fx + fx - fx + h + e + u, h + h -> h + h + e + 9, (3) 
fi + h + e - h + fx + h + H + e - h + h + (4) 

and the modified bremsstrahlung process 

f + f^f + f + v + 9. (5) 

The modified processes are characterized by a spectator baryon that guarantees energy and momentum 
conservation in baryonic matter. In quark matter these processes are subdominant due to the extra 
phase space required by the spectator quarks. Indeed each extra fermion in the initial and final state 
introduces a small factor T/Ep <C 1, where Ep is the Fermi energy. The general arguments above apply 
to the reactions in quark matter featuring strange quarks and in the hypernuclear matter, where the 
kinematical constraints are less restrictive than in purely nucleonic matter [61 [7] . 

Due to the attractive component of the strong interaction nucleons and quarks form Cooper pairs 
at sufficiently low temperatures (for up-to-date reviews on nuclear and quark superconductivity see 
Refs. [U [9]). The formation of Cooper condensates lifts the constraint on the neutral current one-body 
processes in nucleonic |10|, [TT] and quark matter [12j , thus leading to the reaction 

{//}-► / + / + "/ + */, / + /-{//} + + (6) 

where {//} refers to a Cooper pair, / + / to two quasiparticle excitations. These processes - termed 
Cooper pair breaking and formation (CPBF) reactions [13] - are efficient in the temperature domain 
T* <T < T c , where T c is the critical temperature of superfluid phase transition and T* ~ 0.2 T c ; they are 
suppressed asymptotically at low temperatures as exp(— A(0)/T), where A(0) is the zero-temperature 
pairing gap. The temperature domain above matches firmly with characteristic temperatures in the 
neutrino cooling era (T c ~ MeV for nucleonic matter). Thus, the CPBF processes are an important 
ingredient of the cooling of at least the nucleonic matter. The case of quark matter is less clear: the critical 
temperature of pairing of quarks in the dominant pairing channels could be as large as 50 MeV; however 
smaller, ~ keV, gaps were predicted for some combinations of quantum numbers, and the associated 
critical temperatures lie within the relevant temperature range [9]. The neutral current processes © 
induced by the superfluidity have their charged current counterparts j!4^ I15j. While the former vanish, 
when the temperature approaches the critical temperature of superfluid phase transition, the emissivity 
of the latter process approaches the value of the corresponding Urea process. 

2 Quantum kinetics of neutrinos in matter 

Among the methods that are used to compute the rates of neutrino production in dense matter those that 
use the language of many-body theory are particularly suited, as the whole approach can be organized 
in a systematic way, that is consistent with the treatment of related problems of the equation of state, 
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specific heat of matter, pairing fields, etc. In particular, the formulations based on the real-time Green's 
functions (RTG) technique allow for treatments of non-equilibrium processes, including situations far from 
equilibrium. The RTG technique was applied to compute the neutrino emissivities for several reactions in 
nucleonic matter by Voskresensky and Senatorov [H] . In their approach the rates are computed from the 
5-matrix with the help of the optical theorem. Alternatively, the neutrino emission rates can be derived 
directly from a quantum kinetic equation for neutrinos, whereby the collision integrals are expressed in 
terms of neutrino self-energies [HE]. Below, the latter method will be illustrated on a few examples. 



2.1 Transport equations for neutrinos 

We wish to write down a transport equation for neutrinos in a general form involving only Green's 
functions and self-energies. One way of doing this is to start with the Dyson equation for neutrinos 
written on a real time contour. At the first order in the gradient expansion, and upon taking the 
quasiparticle limit in neutrino propagators, one finds [UJ 



{Re S~\q, x), 5 > ' < (g, x)} pB = S>><(q, x)0>'<( 9 , x) + 0>><(q, x)S > ' < (q, x) 



(7) 



where S > ' < (q,x) and Q > ' < (q,x) are the neutrino propagators and self-energies, q = (qo,q) and x = (t,x) 
are the four-momentum and space-time coordinates, {...}p.b. is the four-dimensional Poisson bracket; 
the symbols >, < refer to the positioning of the time arguments of the two-point functions S and O 
on the real-time contour; KeS~ 1 (q, x) is the inverse of the retarded Green's function. The 1. h. side 
of Eq. (J7J) corresponds to the drift term of the Boltzmann equation (hereafter BE), while the r. h. is 
the collision integral. The on-mass-shell neutrino propagator is related to the single-time distribution 
functions (Wigner functions) of neutrinos and anti- neutrinos, f v (q,x) and fp(q,x), via the ansatz 



So(q,x) 



ITT i 



5 {q - uj u (q)) f u (q, x) - 5 (q + uj v (q)) (1 - fp(-q, x)) 



(8) 



where Lo u (q) = \q\ is the on-mass-shell neutrino/anti- neutrino energy. Note that the ansatz includes 
simultaneously the neutrino particle states oc f u (q,x) and anti-neutrino hole states oc 1 — fp(—q,x). 

Upon applying the trace operation (in the space of Dirac matrices) on both sides of the transport 
equation ([7]) and integrating out the off-shell energies on the 1. h. side, one obtains a single time BE for 
neutrinos 



d t + d g u v (q)d x f„(q,x)= I ^-Tr [n < (q,x)S^ (q,x) - n> (q,x)S^(q,: 



0) 



a similar equation follows for the anti-neutrinos if one integrates in Eq. ([7]) over the range [— oo,0]. 



2.2 Collision integrals 

Leading order contributions to the neutrino radiation rates arise from second order Born diagrams when 
the neutrino self-energies f2 >,<; (q) are expanded with respect to the weak coupling constant. The diagrams 
contributing to the charged and neutral current processes are shown in Fig.[TJ The corresponding neutrino 
self-energies are given by 

-Xl>><fa,z) = J ^ (2v)H\ qi -Q2- q)ir^S<(q 2j x)irl\U> k <(q,x), (10) 

where n^ < ((7) refer to the polarization tensors, Tq is the weak interaction vertex to be specified below. 
The central problem of the theory is to compute the polarization tensors of nucleonic or quark matter. 
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Figure 1: The neutrino self-energy in the case of charged (left) and neutral (right) current interaction 
mediated by W + and Z° gauge bosons, respectively; the baryon/quark propagators are labeled as fi, 
i = 1,2; that of neutrino, anti-neutrino and electrons as u, v and e~ . Note that the v and v propagators 
are shown for illustration and should not be included in the evaluation of the diagram. 



To obtain the emissivity through, e. g., a charged current process, we compute from the BE the 
change in the energy per unit volume and time due to the change in the anti-neutrino distribution 

5{uo e +u: u - q Q )u v {q 2 )g B {q ) [1 - f B (u p )] [1 - f e {u e )] A^(qi, q 2 ) Im II* (q), (11) 

where G = Gcos 9c, G is the weak coupling constant, 9c is the Cabibbo angle (cos 9c = 0.973) and 
(<7i, q 2 ) = Tr [7^(1 — 7 5 ) )^i7 < '(l — 7 5 ) gfe] • The symbol Im refers to the imaginary part of the po- 
larization tensor's resolvent. Here we used the relation IT^(g) = 11^ (—q) = 2igB{qo)ImIl^(q), where 
9b{<1o) is the Bose distribution function and H^(q) is the retarded component of the polarization tensor. 
In equilibrium, /p/ e (wp/ e ) reduce to Fermi-distribution functions for anti-neutrinos and electrons. Since 
the anti-neutrinos leave the star without interactions, there is no thermal population of anti-neutrinos, 
i. e. fu{uJu) <ti 1 and can be neglected. The neutrino emissivity for the case of neutral current processes 
can be obtained in a similar way [H |T7j . 



3 Polarization tensors of dense matter: Examples 

It is instructive to study the polarization tensors describing charged and neutral current processes first 
at the single loop level. Descriptions that are consistent with the conservation laws and Ward identi- 
ties require vertex corrections to the one-loop results, which we shall address later on. We shall now 
switch to the equilibrium finite temperature techniques of Matsubara Green's functions thus treating the 
nucleonic/quark matter in thermal equilibrium. 



3.1 Direct Urea process in baryonic matter 

Since the temperature of dense matter core during the neutrino cooling era is well below the critical 
temperatures of pairing in baryonic and quark matter, pairing correlations should be included in the 
computation of polarization tensors. The non-relativistic Matsubara propagators that incorporate the 
pairing correlations are given by 

/ u 2 v 2 \ 

Gt'(iPn,p) = 8 aa ,6 TT , . v _ + . * = 5 a(7l 5 TT ,G M (ip n , P ), (12) 

F aa '(iPn,P) = -i(T y 6 TT 'U p V p ( - — — ), (13) 

\iPn - Bp iPn +£pj 
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where p n = (2n + 1)tvT is the fermionic Matsubara frequency, a and r refer to spin and isospin, a y 
is the y-component of the Pauli-matrix, u 2 = (1/2)(1 + i p /e p ) and v 2 = 1 — u 2 are the Bogolyubov 

amplitudes and e p = W 1 Z 2 + A| is the quasiparticle spectrum, where £ p = p 2 /2m* — p is the spectrum 

in the unpaired state, with m* and p being the effective mass and chemical potential. Here A p is the 
anomalous self-energy (gap function). The propagators above are written for the case of S- wave neutron 
or proton proton pairing in isospin-1, spin-0 state. Note that at high densities the neutron fluid is paired 
in a P wave (this is not the case for protons because of their low abundance). Consider now the case of 
direct Urea process involving nucleons (Fig. [IJ left diagram). The Matsubara polarization tensor is then 
given by 

mfKiq, q) = J2[ 4^ Tt [r^L 1/2 (ip,p)r^L_ 1/2 (ip + i q , P + g)] , (u) 

ip 

where the charged current weak interaction vertices are = 7^(1 — £L47 5 ), with qa = 1.26 being the 
axial coupling constant. Upon performing the Matsubara sums and analytical continuation (iq — > iv + i5) 
we obtain the retarded polarization tensor 

tR <A„ „\ — O / & ) I "p^p+q u p u p 




n$ /A (iq, q ) = 2 / ^- * ™ lf(ep) - f{e p+q )} 

1 J X^) 6 y yu + e p -s p+q +ib lo - e v + £p +q + id J 

' p ^ \LU - e p - £ p+q + lb UJ + £ p + e p+q + id J J 

where the vector/axial- vector polarization tensors are the components proportional to 1 and g\, 

respectively. The first two terms in Eq. (|15p correspond to excitations of a particle-hole pair while the 
last two to excitation of particle-particle and hole-hole pairs. The last term does not contribute to the 
neutrino radiation rate {to > 0). We identify the first two terms as the scattering (SC) terms, while 
the third term as the pair-braking {PB) term. Upon evaluating the phase space integrals, the neutrino 
emissivity is written as ep = eQ rca J, where 

^o Urca = (i + 3 g i) 3G2< 2 7f FeT6 , J = -lfj ygB (y)[i sc (y) + i PB (y)]fdzz*Uz-y), (16) 

where pFe is the Fermi-momentum of the electrons and y = j3u; the integrals I sc {y) and I PB (y) are 
given in Refs. [14J. In the unpaired state (u p — ► 1 and v p — > 0) only the scattering contribution survives; 
upon integrating we obtain 



I SC (y) = In 



1 + exp [-PZ] 



1 + exp [-/?(£ + u) 



(17) 



ei 



where ^ = p 2 /2m* — /x p and p = (m*/q)(u — [i p + \i n — q 2 /2m*); here the momentum transfer q = pp 
fj, n and Hp are the chemical potentials of neutrons and protons, and we assumed for simplicity that their 
effective masses are equal. In the zero temperature limit I sc (y) = y0(—/3£), the integrals in Eq. (|16() 
can be performed analytically and one recovers the zero-temperature result of Lattimer et al. [2]. The 
zero temperature ^-function can be rewritten as 0(pp e + pp p — pFn) [2] which tells us that the "triangle 
inequality" PFe+PFp > PFn must be obeyed by the Fermi- momenta of the particles for the Urea process 
to operate. 

While the one- loop approximation provides a useful starting point, the complete treatment of the 
problem when the particle-hole interaction is not small, i. e. can not be treated as a perturbation, 
requires summation of infinite series of particle-hole loops. This is certainly the case in nuclear matter, 
where the Landau parameters are 0(1). Figure [2] shows the four distinct diagrams in the case where the 
loops are summed up to all orders. Next subsection shows how to improve on the one-loop result using 
as an example the neutral current interactions. 
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Figure 2: The sum of polarization tensors that contribute to the neutrino emission rate. The single 
arrowed solid line corresponds to the normal propagator (|12p . while the double arrowed line to the 
anomalous propagator f)13[) . The contributions form H^ b \q), H^ c \q), and U^ d \q) diagrams vanish at 
one-loop for the Urea process; the U( d \q) dia gram is finite for charge neutral interactions at one- loop. 



3.2 Neutral current neutrino pair-bremsstrahlung 

The polarization tensor describing the neutral current interactions in baryonic matter is given by 

ilV (9) = / Tr [ T - G ^ P) T MiP + iq,P + q)+ T u F(ip, p)T^ {ip + iq,p + q)] , (18) 

ip 

where the neutral current vertices are = 7^(cy — C/17 5 ). Performing the Matsubara sums we obtain 
for the vector and axial-vector contributions in this case 

n^(g) = Z [f(e p ) - f(e k )] (-— ^ . \ ) 

^ \tq + e p -e k iq - e p + e k ) 



CV Dn 



where k = p + q, A T = u 2 u\ =p h, = v 2 v\ =F h, C T = u\v 2 ± h, D T = u 2 v v\ ± h, h = u p u k v p v k . 
The first line in Eq. (|19p corresponds to the process of scattering where a quasiparticle is promoted out 
of the condensate into an excited state, or inversely, an excitation merges with the condensate. The 
corresponding piece of the response function ImU v / A (q) vanishes for small momentum transfers. The 
second line in Eq. (119p describes the process of pair-breaking and recombination, i. e., excitation of pairs 
of quasiparticles out of the condensate, and inversely, restoration of a pair within the condensate. Since 
we are interested in the emission process we shall keep only the terms that do not vanish for u > 0; then, 
the pair-braking contribution is given by the term oc C±. This contribution to the polarization tensor 
can be evaluated analytically in the limit q —* and the case A 7^ A(p) and is given by 



(20) 



ImnVfe) = -2^p F H^fQ 2 ^) 7 =^=e(.-2A), 

lmU A (q) ~ + O f^fV (21) 

where v(jPf) = rn*pp/2-K 2 is the density of states (h = 1) and 9 is the Heaviside step function; the explicit 
form of the O (u^/c 2 ) contribution to the axial current response is given by Flowers et al in Ref. |10j . 

Upon substituting Eq. ([20]) in the neutral current analog of Eq. (fTT|) and carrying out the phase-space 
integrals we obtain the emissivity per neutrino flavor [10] 

' - °~ <T Hpf) T 7 h(C) = 6^ rcms - h(0, (22) 



vv 240vr 3 
where C = 2A(T)/T and 

h(C) = C 



J d4> (cosh 4>f f f| cosh . (23) 
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Ta ^\Ai = ^C.^-^^^^ +^^*^^^\/v> +^^^^^xx\» +^^^^^^>^ 

Figure 3: Coupled integral equations for the effective weak vertices in superfluid baryonic matter. The 
"normal" T± vertex (full triangle) and two "anomalous" vertices T2 (hatched) and T3 (shaded triangle) 
are shown explicitly, the fourth vertex (empty triangle) is obtained by interchanging the particle and hole 
lines in the first line. The anomalous vertices vanish in the normal state. 

Note that the rate (|22|) scales as £ and, consequently, it is sensitive to the magnitude of the pairing gap. 
Because of the substantial density dependence of the gap, the emissivity ([22]) varies strongly across the 
stellar interior. 

In nuclear and neutron matter problem the particle-hole interactions are not small and cannot be 
treated in the perturbation theory. The resummation of particle-hole diagrams in a superfluid leads to 
coupled integral equations shown in Fig. [3l In the non-relativistic limit the driving terms in the vector 
and axial-vector channels correspond to the scalar and spinor perturbations, i. e. the bare vector and 
axial- vector vertices are T v = 1 and T A = er. The topologically non-equivalent polarization tensors and 
the associated vertices are shown in Figs. [2] and El 

Including the vertex corrections modifies the one-loop result to (Sedrakian et al. in ref. |10| ) 

evp = e (h + h), (24) 

f'OO ^ 

h = r^7 du u 6 g{u)J2[^i^) ReTiM], (25) 

2 f'OO 3 

h = r^f duJUJ & g{oj) V[ReIT(u;) ImT^o;)] , (26) 



where 



where IIi = II^ - II W , U 2 = 11^ and II3 = 11^ . For T -> T c the rates vanish, consistent with 
the observation that the pair bremsstrahlung is absent in normal matter for on-shell (non-interacting) 
baryons. At small T < 0.3T C the rates are suppressed exponentially as exp(— A/T). 



3.3 Direct Urea process in color superconducting matter 

Although quark matter in compact stars and its superfluidity were suggested more than three decades ago, 
these topics have received much attention in recent years after the models, which were designed to describe 
the chiral phase transition in dense matter, were applied to the problem of quark superconductivity (the 
current state of the art is reflected in the reviews [9]). At moderate densities relevant to compact stars 
the quark matter is in the non-perturbative regime and one has to rely on effective models that capture 
(at least some) features predicted by the QCD (chiral symmetry breaking, confinement, etc.) The ground 
state of superconducting quark matter under /^-equilibrium is not known; one significant problem is that 
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Figure 4: One-loop ^-polarization tensor: (a) normal piece, which represents the sole contribution 
above T c ; and (b) "anomalous" piece, which is proportional to A 2 (T) and vanishes at T c . 



under the stress caused by /3-equilibrium and/or the strange quark mass, the Fermi-surfaces of up (u) 
and down (d) quarks are shifted apart, and the resulting pairing patterns differ from the ones predicted 
by the BCS theory. Color and flavor degrees of freedom are responsible for the multitude of possible 
pairing patterns. 

In contrast to the case of nucleonic matter the Urea process is permitted in interacting quark matter 
for any asymmetry between u and d quarks. The emissivity to first order in the strong coupling constant 
a c is given by [3] 

e o = TT7 G 2 F ol c fi d fj, u fi e T 6 , (27) 
315 

where fa, i = d, u, e are the chemical potentials of down and up quarks and electrons. If non- 
superconducting quark matter is present in the core of a neutron star, the star will cool very rapidly to 
temperatures well below the observational threshold. 

What are the effects of superfluidity on the cooling rate of such a star? Let us consider a specific 
model [15] . where the pairing is in the so-called 2SC phase, i. e. quark pairing is characterized by the 
order parameter A oc (ip T (x)C~/5T2\2ip( x )} > where T2 is the Pauli matrix in the isospin state, A2 is the 
Gell-Mann matrix in the color space, C = ij 2 ^/ is the matrix of charge conjugation. The minimal 
effective Lagrangian describing the pairing is given by 

£ cfr = i>{x){i^d^{x) + G^C^Xa^x^^C^Xa^x)), (28) 

where G\ is the attractive pairing interaction. The Lagrangian is minimal in the sense that apart from 
the pairing interaction it includes only the kinetic term for massless quarks; other channels of interaction, 
such as the repulsive components which would lead to the renormalizations of the single-particle spectra 
of quarks are omitted. The normal and anomalous propagators of quarks of flavor / are 

A + (p) A + (p) 

s f=u , d = ^ (PD + W _^ »~ Hp) = -» eaMe/g A (pD + w _ eg 75 c, (29) 

where e p = \/ (p — fi) 2 + A 2 , 5fi = (fid — Mu)/2 and \i = (fid + ^ u )/2, A + is the projector to the positive 
energy state. The emissivity at one-loop can be obtained by evaluating the sum of diagrams in Fig. [U 
which leads to the polarization tensor 

n,x(q) = -i J ■0^TT[(T^) ll S(p)(T + ) x S(p + q ) + (T^ l ,F(p)(T + ) x F(p + q )}, (30) 

where r-t(g) = G"f^(l — 75) ® t± and t± = (t\ ± T2)/2 are flavor-raising and lowering operators. The 
result for the emissivity depends on whether the parameter £ = A/5fA is larger or smaller than unity [15] . 
For £ > 1 all the particle modes are "gapped", therefore, as the temperature is lowered, the emissivity is 
suppressed (for asymptotically low temperatures exponentially). When £ < 1 there are gapless modes in 
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Figure 5: Temperature dependence of the neutrino emissivity normalized to its value at the critical 
temperature for a range of values of ( = A(0)/5fi; the upper panel corresponds to the gapped regime 
(C > 1) while the lower panel to the gapless regime (£ < 1). The dashed and dashed-dotted lines are the 
normal and anomalous contributions; the solid line is their sum. 

the quasiparticle spectrum; this implies that the neutrino production is not affected by color supercon- 
ductivity for these modes. As a result, the superconducting quark matter cools at a rate comparable to 
the unpaired matter. Fig. [5] illustrates these two distinct cases. 

The neutrino emissivity of color-superconducting matter has been studied for alternative realizations 
of the ground state matter: one such realization is the spin-1 color superconductivity |16| . Since the 
condensate in this case breaks the rotational symmetry the neutrino emission turns out to be anisotropic 
for some choices of the order parameter [18]. Another realization is the crystalline color superconduc- 
tivity, which is characterized by spatially modulated gap parameter. The neutrino radiation rates from 
crystalline color superconducting matter and cooling of compact stars featuring such a phase is discussed 
in Refs. p2]. 

4 Suggested exercises 

1. Derive the emissivity of the direct Urea process n — > p + e + v in unpaired matter by using the 
Fermi Golden rule and following the similar derivation of the emissivity of the modified Urea process 
n + n — >n + p + e + v'm. Ref. [20] . Rederive the emissivity in unpaired matter by setting u p = 1 
and v p = in Eq. (JTSD and A^ImlT^ ~ 8u e uj u (lmll v + 3g\lmIl A ) in Eq. (TTT]) . 

2. Derive the emissivity of unpaired quark matter featuring u and d quarks through the direct Urea 
process d — ► u + e + vhy using the Fermi Gold rule in the case where the quarks interact to leading 
order in a s [3]. Repeat the calculation by starting from the polarization tensor (I30|) with A = (see 
Ref. [15J for the case A / 0). Next assume that quarks are non- interacting but massive. Repeat 
the calculations in this case and compare to the result of Ref. [3] . 
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